Spin relaxation in n-type ZnO quantum wells 
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We perform an investigation on the spin relaxation for n-type ZnO (0001) quantum wells by nu- 
merically solving the kinetic spin Bloch equations with all the relevant scattering explicitly included. 
We show the temperature and electron density dependence of the spin relaxation time under various 
conditions such as impurity density, well width, and external electric field. We find a peak in the 
temperature dependence of the spin relaxation time at low impurity density. This peak can survive 
even at 100 K, much higher than the prediction and measurement value in GaAs. There also exhibits 
a peak in the electron density dependence at low temperature. These two peaks originate from the 
nonmonotonic temperature and electron density dependence of the Coulomb scattering. The spin 
relaxation time can reach the order of nanosecond at low temperature and high impurity density. 



PACS numbers: 72.25.Rb, 73.21.Fg, 71.10.-w 
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I. INTRODUCTION 

Much attention has been devoted to the spin degree 
of the freedom of carriers in the zinc oxide (ZnO) with 
wurtzite structure in the last few years, ^ partly because 
of the very long spin relaxation time (SRT)^i^ii and 
the prediction that ZnO can become ferromagnetic with 
a Curie temperature above room temperature if doped 
with manganese^^ However, only few works investigate 
on the spin dynamics properties of ZnO: Experimentally, 
Ghosh et al^ investigated the electron spin properties in 
n-type bulk ZnO and discovered the electron spin relax- 
ation time varying from 20 ns to 190 ps when the temper- 
ature increased from 10 to 280 K. Liu et al.^ measured 
the SRT in colloidal 7i-type ZnO quantum dots could be 
as long as 25 ns at room temperature by electron param- 
agnetic resonance spectroscopy. Theoretically, Harmon 
et al. calculated the SRT in bulk material in the frame- 
work of a single-particle modelj^ However, a theoretical 
investigation on the SRT in quantum wells (QWs) is still 
rare, which is our task in the present paper. 

The spin relaxation can be induced by the following 
mechanisms: (i) D'yakonov-Perer(DP) spin relaxation 
mechanism,^ which is revealed to be the dominant mech- 
anism in n-doped semiconductors.^'^ For n-doped ZnO 
QW, the spin-orbit coupling (SOC) is at least one order 
of magnitude weaker compared with that of the well stud- 
ied semiconductor GaA a^°i^^ , a simple estimation shows 
that the criterion of the strong scattering^ is always sat- 
isfied even the momentum scattering is also weaker than 
that in GaAs due to the larger effective mass m* . There- 
fore the DP mechanism here can be described by the mo- 
tional narrowing picture^^ qualitatively and the induced 
SRT is 
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SOC. As pointed out first by Wu et alM-^^ and then by 
Glazov and Ivchenkofi^ the electron-electron scattering 
has important contribution to the spin relaxation pro- 



cess. Therefore r„ 



as T,: 



used in Eq. 
1 



ID) should be revised 



to include the electron- 
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with Tp standing for the momentum scattering time and 
(il^) for the inhomogeneous broadening induced by the 



electron momentum relaxation time Tp^.- 
(ii) Elliott- Yafet mechanism. The revised criterion in 
Eq. (22) of Refs. ;8J gives 9 « 340 eV, which means the 
Elliott- Yafet mechanism is negligible compared with the 
DP mechanism^ due to the small spin split off energy, the 
large band gap, and the large m*. (iii) Bir-Aronov-Pikus 
mechanism)^ which is always unimportant in n-doped 
semiconductor!^ Therefore, we only investigate the SRT 
induced by the DP mechanism for n-type ZnO QWs in 
the following. 

In this paper, we quantitatively calculate the SRT for 
n-type ZnO QWs by using the fully microscopic spin ki- 
netic Bloch equation (KSBE) approach, which has been 
successfully used in investigating the spin relaxation in 
Q-^99ii3ii4ji5^ and in bulk semiconductors,2i2^ With all 
the relevant scattering included, the influence of temper- 
ature, electron density, impurity density, well width and 
electric field on the SRT are studied detailedly. The tem- 
perature and density dependence of the SRT is shown to 
be nomonotonic, and we find that the SRT increases with 
the electric field monotonically. 

This paper is organized as follows: In Sec. 11 we de- 
scribe our model and the KSBEs. Our numerical results 
are presented in Sec. 111. We conclude in Sec. IV. 



II. MODEL AND KSBES 

We start our investigation from a n-doped ZnO QW 
of well width a grown in (0001) direction, considered to 
be z axis. Due to the confinement of QW, the momen- 
tum states along 2:-axis is quantized by subband index 
n. With the momentum vector k = [k. 



ky) and the 



2 



spin index a, the electron Hamiltonian can be written 



as H„ 



nai (72 k 



Here Sncr-ia^iz 



Erik 



2m» 



= e,ik(5criCT2 + [hj?,(k) + (h_D)„(k)]-crCTjcr2 with 
is the energy spectrum, R = (a;,?;) 



is the position, the effective magnetic field given by the 
Rashba due to the intrinsic wurtzite structure inversion 
asymmetry and Dresselhaus SOC can be written asJ^ 



hfl(k) = ae{ky, -fcx,0) , 



(hz,)„(k) ^ 7e{b{kl), 



fc^)(fcj^. 



(2) 



with «£, 7e and b standing for the SOC coefficients. 
{kl)n = ^ ^ii" is the subband energy in a hard- 
wall confinement potential. The scattering Hamilto- 
nian Hj includes all the scatterings, such as electron- 
nomagnctic impurity scattering, electron-phonon scatter- 
ing, and electron-electron scattering. 

We construct the KSBEs in the coUinear statistics 
by using the non-equilibrium Green function method as 
foUows-^^i^s^ 



dtpw - eE • VkPk = <9f/9k|coh + dtPk\ 



(3) 



The density matrix pk for momentum k is a matrix with 
matrix elements [pk]nia-i-n2iT2 which include all the coher- 
ence between different subbands and different spins. The 
second terms on the left-hand side of the kinetic equa- 
tions describe the electric field E driven effect. i9tPk|coh 
is the coherent term. i9tPk|scat denotes the scattering, 
including the electron-impurity, the electron-phonon, as 
well as the electron-electron scattering. The expressions 
for these terms are given in Appendix A. 

Before we give our numerical results, the qualita- 
tive analysis of the DP mechanism due to the electron- 
electron scattering can be made at strong scattering 
limit. The perturbation theory shows the effective 
electron-electron momentum scattering time Tp*^ in de- 
generate and nondegenerate limits satisfies^^ 
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for T < Tf 
for T:$>Tf 



(4) 



which has nonmonotonic temperature T and electron 
density Ne dependence as the electron gas undergoes the 
transition from the degenerate case to the nondegenerate 
case at the Fermi temperature Tp. As the SOC in ZnO 
QWs mainly depends on k linearly, the inhomogeneous 
broadening is given by 



Ne, for T^Tf , 
T, for T:^Tf . 



(5) 



Then by Eq. ([T]), the electron-electron scattering con- 
tributes to the SRT T atM 



T^N-^, for T^Tf , 
T~^Ne, for T:^Tf . 



(6) 



The SRT is expected to reach a minimum in T depen- 
dence or a maximum in iVg dependence, and the location 
of the extreme points satisfie T ^ Tf- 



III. NUMERICAL RESULTS 

We numerically solve the KSBEs for the spin density 
matrix p, from which we obtain the time evolution of the 
spin polarization along z-directions^^ 



Pzit) = $^{[pk]nT;nTW - [Pk]ni;niW}/iVe 



(7) 



where Ne is the total electron density. The SRT t is 
extracted from the exponential decay of the envelope of 
Pz (t) . The initial condition at t = is taken to be 



(8) 



Here /T,/.(k) = [1 + e(^''-M)/(feBT)]-i gj^gg Fermi- 
Dirac distribution. The spin dependent chemical po- 
tential Per is chosen to satisfy P(0) = 2.5%. The elec- 
tron density and the quantum well width are taken as 
A^e = 4 X fO^^/cm'^ and a = 10 nm respectively unless 
otherwise specified. All used parameters are listed in Ta- 
ble m In the calculation, only the lowest two subbands 
are taken into account. 



TABLE I: Material parameters used in the calculation (from 
Ref . [2^ unless otherwise specified) . 



/^oo 3.T 

nie/mo 0.25 

V3 (km/s) 6.09 

va (km/s) 2.79 

7e (eVA^) 0.33" 

ei5 (V/m) -0.35 x lO'^ 



Ko 7.8 
vi (km/s) 6.08 
«4 (km/s) 2.73 
b 3.91" 

(meVA) LI*" 
633 (V/m) 1.56 X 10^ 



Ref. [3; 



Ref. hJ 



A. Temperature dependence 

We now study the temperature dependence of the SRT 
presented in Fig. [T] for different impurity densities. The 
results are similar to that in GaAs QWs^'^^ and can be 
understood as follows: (i) The SRT always increases with 
the impurity density A'^^. It is because the system is 
in the strong scattering regime as stated above due to 
the weak SOC, thus the spin relaxation can be explained 
by the motional narrowing picture qualitatively^iii^, and 
the additional scattering leads to longer SRT. (ii) The 
electron-phonon scattering is shown to be negligible over 
the whole temperature regime by comparing the tem- 
perature dependence of SRT with (solid curve with A) 
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FIG. 1: SRT r vs. temperature T at difTerent impurity den- 
sities. The dashed curve is obtained from the calculation of 
excluding the electron-phonon scattering. 



FIG. 2: SRT vs. the electron density with different impurity 
densities and temperatures, (a) T = 20 K; (b) T = 300 K. A: 
N,/Ne = 0; ■: Ni/Ne = O.I; N,/N, = 1. 



or without (dashed curve with A) the electron-phonon 
scattering for the impurity free case in the same fig- 
ure, (iii) The SRT presents a peak at very low Ni. In 
these cases, the electron-electron scattering is the domi- 
nant scattering, therefore, as shown in Eq. ([6]), the SRT 
shows a maximum and the transition temperature Tp 
for iVe = 4 X 10^^ /cm^ is about 44 K, which is agree 
with our numerical results, (iv) When the impurity den- 
sity is high enough such as Ni = Ne, the SRT decreases 
monotonically with T. In this case the total scattering 
is mainly determined by the impurity scattering, which 
depends weakly on the temperature. However, the inho- 
mogeneous broadening from the DP term increases with 
the temperature, and results in shorter SRT. 

For GaAs QWs, the temperature peak of the SRT can 
only be observed at low electron density (i.e. low transi- 
tion temperature) and low impurity densitySi2£j because 
the electron-phonon scattering becomes strong enough to 
destroy the nonmonotonic T dependence of the scattering 
time induced by the electron-electron scattering. Such 
case can be avoid in ZnO QWs, in which the electron- 
phonon scattering is always pretty weak due to the large 
optical phonon energies (^ 800 K). Thus the temperature 
peak can be found even for high electron density samples. 
One can easily find from the chained curve in Fig. [1] 
which is calculated with parameters Ne = 10^^ /cm^ and 
Ni = 0, that the peak moves to T ~ lOOK. Therefore, 
the high mobility ZnO QW is a good system for studying 
the electron-electron scattering. 



B. Doping and well width dependence 

Then we investigate the density dependence of the 
SRT at different temperatures and impurity densities. 
In Fig. 2 (a) we plot the SRT as a function of the 
electron density with T — 20 K. One can see that for 
the low impurity density case, the SRT reaches a max- 
imum at A^e ~ 2 X 10^-'^ cm~^, which has been pointed 
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FIG. 3: SRT vs. the temperature at different quantum well 
widths. 



out in n-type bulk III-V semiconductors^i^^. It origi- 
nates from the transition from the nondegenerate elec- 
tron gas to the degenerate electron gas and can be well 
explained by Eq. ([6]). Our calculation gives the transi- 
tion density of 2 x 10^^ cm~^, corresponding Tp ^ 22 K, 
close to the lattice temperature of 20 K. For the case 
of Ni — Ne, the electron-impurity scattering time has 
the same Ne dependence as that for electron-electron 
scattering: in the nondegenerate regime, oc Ni{Ug), 

in which (U^) changes little; in the degenerate regime, 
~ NiU^ , oc Ne/kp (X Ne^. Consequently, the peak 

still exists and is almost at the same position. In com- 
parison, the SRT as a function of Ne with T = 300 K is 
plotted in Fig. 2 (b). In this case, one finds that the SRT 
increases monotonically with A^e. This could be easily 
understood for that Tp ^ T is satisfied and it is in the 
nondegenerate regime, in which the SRT increases with 
density as discussed above. 

We further show the effect of quantum well width on 
the spin relaxation. In Fig. 3 the SRTs versus temper- 
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FIG. 4: SRT vs. the electric field at different temperature 
and impurity densities. A: N^/N^ = 0; ■: Nt/Ne = 0.1; •: 
N^/N, = 1. 



temperature. Both these two peaks originate from the 
different temperature and electron density dependence 
of Tp"^ in degenerate and non-degenerate case. Compared 
with the same effect in III-V semiconductor this 
peak position can occur at the temperature as high as 
100 K and is easier to observe in experiments due to 
the weak electron-phonon scattering. When the impu- 
rity density is high, the peak in the temperature de- 
pendence disappears and the SRT decreases with tem- 
perature monotonously. Moreover, the peak in the elec- 
tron density dependence moves to larger electron density 
which is beyond the scope of our interest when the tem- 
perature is high. We also investigate the hot-electron 
effect and show that the SRT always increases with the 
electric field. It is also shown that the SRT reaches the 
order of nonosecond at low temperature and high impu- 
rity density. 



ature at well widths a = 10 nm and 20 nm are plotted 
respectively. Both the SOC and the scattering^^'^^ de- 
pend on the quantum well width. However, comparing 
to the weak well width dependence of the scattering, the 
fast decrease of (/c^) in the DP term with a dominates 
and so the SRT increases with well width. 



C. Electric field dependence 

Then we investigate the electric field dependence of the 
SRT at different temperatures and impurity densities. In 
Fig. 4 we plot the SRT as a function of the electric field for 
different T. The electric field is applied along the x axis. 
One can see that the SRT increases monotonically for 
both low temperature and high temperature cases. Ac- 
cording to the previous investigation^- , the electric field 
will enhance both the momentum scattering due to the 
hot-electron effect, and the inhomogencous broadening 
due to the drift of the electron distribution to larger k 
states. These two effects are competing effects for the 
SRT: the former tends to enhance the SRT while the 
later tends to suppress itii^ For SOC with hnear k de- 
pendence, the hot-electron effect dominates,^ thus the 
SRT always increases with E. 



IV. CONCLUSION 

In conclusion, we have investigated the spin relax- 
ation for n-type ZnO (0001) QWs by numerically solving 
the KSBSs with all the relevant scattering explicitly in- 
cluded. It is shown that the electron-phonon scattering is 
pretty weak in ZnO QWs, while the Coulomb scattering 
always plays an important role. Therefore the ZnO QW 
is a good carrier to study the electron-electron scattering. 
We find there exists a peak of SRT both in the tempera- 
ture dependence for a given electron density at low impu- 
rity density and in the electron density dependence at low 
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APPENDIX A: EXPRESSIONS FOR KINETIC 
BLOCK EQUATIONS 

Here we write the expressions for the coherent terms 
and the scattering terms in the kinetic Bloch equations. 
The coherent terms in Eq. ([3]) can be written as 



dtPk 



coh 



Q 



(Al) 



where [A, B] = AB — BA denotes the commutator. 



ni(7i ;r 



'^^ni(Ticr2k<5niri2 ■ The Coulomb Hartree- 



Fock term, which is always negligible for small spin 
polarization,— i2i is also included. The form factor Iq^ 
is also a matrix with matrix elements 



[I 



qJni(Ti;n2cr2 

1 



= W(7(5cri,(T2 [s'"' COS7r(71i — ^2) — 1] 

1 



(A2) 



The statically screened Coulomb potential in 
the random-phase approximation (RPA) reads^^ 
Vq = with the dielectric function e(q) = 

1 _ V l)r^\\T 1 p /^(g.il,k+q)-/,^(£n2,k) , 



Q = {ixilyilz)- The bare Coulomb potential is 
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The scattering terms in Eq. ^ can be written as limit with 
StPklscat = dtPklim + ^tpklph + ^tPkloG in the Markovian 



dtPk 
dtPk 

dtpv 



ph 



£rii,k-q £n2,k 

)h. [(1 - Pk-q)TnJ-q,Tn2Pk - Pk-qTnJ-q.Tn^il - Pk)] + h.C. , 

Q.ni ,712 

E l^'^Q,Ap^9,{<5(eni,k-q - en2,k +t^Q,A)[(^Q.A + 1)(1 - Pk-qjTnJ-q.Tn^Pk 
Q,ni ,712, A 

- A^Q,APk-qr„i/_q^T„2(l - Pk)] + <5(£ni,k-q " £712, k " ^Q,a) [A^Q, A (1 - Pk-cijTnJ^q^Tn^Pk 

- (-^Q.A + l)/9k-qr„i/_q^T„2(l - pk)]} + h.C. , 

TT E E VQVQ'(5(£„i,k-q - en2,k + e>l3,k' - e7l4,k'-q)J'g, 

q,gsg^,k' 711,712,713,714 

X {(1 - Pk-q)T'„i/_g^r„2PkTr[(l - pk')TnsIq'^TniPk' -ql-q'J 

- pk^qTnJ^q^Tn^il - Pk)Tr [pk'T„3 7,/^ T„^ ( 1 - pk'-q)/-gj} + h.C. , (A3) 



in which [T„J„,„' = ^7ii,7i<5„i,7i' and Q' = {qr^,qy,q'J. Ni 
is the density of impurities, and |C/qP is the screened im- 
purity potential. |Mq^a|^ andiVq^A = [ey^Pi^Q,\/kBT) - 
are the matrix elements of the electron- phonon in- 
teraction and the Bose distribution function respectively. 
^Q.A is the phonon energy spectrum. For the electron- 



phonon scattering, we only include electron AC-phonon 
scattering^ for which the explicit expressions can be found 
in Refs. 25 and '26*. The electron-optical phonon scatter- 
ing is ignored for that the optical phonon energies are 
around 70 meV and are out of our temperature rangci^ 
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